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The Gari-Kriimpelmann (GK) models of nucleon electromagnetic form factors, in 
which the /3, and (p vector meson pole contributions evolve at high momentum 
transfer to conform to the predictions of perturbative QCD (pQCD), was recently 
extended to include the width of the p meson by substituting the result of dispersion 
relations for the pole and the addition of the p' (1450) isovector vector meson pole. 
This extended model was shown to produce a good overall fit to all the available 
nucleon electromagnetic form factor (emff) data. Since then new polarization data 
shows that the electric to magnetic ratios Rp and Rn obtained are not consistent 
with the older Gep and Geu data in their range of momentum transfer. The model is 
further extended to include the uj' (1419) isoscalar vector meson pole. It is found that 
while this GKex cannot simultaneously fit the new Rp and the old Geu data, it can 
fit the new Rp and Rn well simultaneously. An excellent fit to all the remaining data 
is obtained when the inconsistent Gep and Geu is omitted. The model predictions 
are extended beyond the data, if needed, to momentum transfer squared, Q^, of 8 
GeVVc^. 

PACS numbers: 13.40.Gp, 21.10.Ft 

I. INTRODUCTION 

A variety of related models of the nucleon emff [|^ were fitted to the complete set of 
data available before September 2001. One group of models included variants of the basic 
GK model of p, a;, and vector meson pole terms with hadronic form factors and a term 
with pQCD behavior which dominates at high [0. Four varieties of hadronic form factor 
parameterization (of which two are used in 0) were compared. In addition to the GK 
type models we considered a group of models (generically designated DR-GK) that use 
the analytic approximation of 0] to the dispersion integral approximation for the p meson 
contribution (similar to that of [§), modified by the four hadronic form factor choices used 
with the GK model, and the addition of the well established p' (1450) pole. Every model 
had an electric and a magnetic coupling parameter for each of the three pole terms, four 
"cut-ofF' masses for the hadronic form-factors and the QCD scale mass scale, Aqcd for the 
logarithmic momentum transfer behavior in pQCD. In addition the effect of a normalization 
parameter was sometimes considered for the dispersion relation behavior of the p meson in 
the DR-GK models. 

When the set of parameters in each of the eight models was fitted to the full set of 
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data available before publication, for Gep, Gmp, Geu, Gmh and the lower values of 
Rp = fipGEp/GMp, three GK and all four DR-GK models attained reasonable (when the 
inconsistency of some low Geu and Gmh data was taken into account), but the extended 
DR-GK models had significantly lower x^. Furthermore Aqcd was reasonable for three of 
the DR-GK models but for only the one of the GK models that had an unreasonably large 
anomalous magnetic coupling Kp. It was concluded that the three DR-GK models were the 
best nucleon emff to use in prediction of nuclear electromagnetic properties. All thee were 
found to be moderately consistent in their predictions up to of 8 GeV^/c^. 

However the part of the above data set from a recent Rp ratio data |^ for 0.5 GeV^/c^ < 
< 3.5 GeV/^/c^, swamped statistically by all the other data, was systematically lower 
than the fitted models (Fig. 5 of contributing disproportionately to x^- This ratio is de- 
termined by an asymmetry measurement in the scattering of polarized electrons on protons. 
Multiplied by the well determined values of Gmp one obtains values for Gep which are not 
subject to the uncertainty inherent in the Rosenbluth separation measurements in which 
Gep is obtained by subtracting the much larger contribution of Gmp from the unpolarized 
cross section. As expected the Gep derived from the measured Rp are consistently below 
those of the older Rosenbluth separation values. 

It is plausible to expect that the old Gep data is responsible for restricting the best fit 
of the models to be substantially above the experimental Rp values. With this in mind the 
particularly high data of was omitted from the fit to the model type DR-GK'(l) of 
and the flexibility of a p meson dispersion integral normalization parameter N was included. 
In this article the original version is designated as GKex(Ol) and when fitted to the smaller 
data set as GKex(Ol-). As seen in Tables | and ^ and Figs. |I| and |^, there is only a small 
change in the fit to Gep and Rp, although the parameters of the fit change substantially. 
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FIG. 1: Gep normalized to Gd, comparing the GKex(Ol) fit [dotted] with the fit GKex(Ol-) [solid] 
obtained when the data of @ is omitted. The other Gep data is from 0, [g, g, |l|], and flit] The 
data symbols are listed in the figure. 



After the publication of ][I| new data [|T3[ extended the measurements of Rp up to = 5.6 
GeV^/c^, exacerbating the discrepancy with the predictions of the best models in |l]]. Very 
recently i?„ = fi'nGEn/GMn has been obtained directly [jl4| by the scattering of polarized 
electrons on deuterium and detecting the polarized recoil neutron at = 0.45, 1.15 and 
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FIG. 2: i?p, the ratio l^pGEp/GMp, comparing the GKex(Ol) fit [dotted] with the fit GKex(Ol-) 
[solid] obtained when the data of [Q is omitted. The data is from [|5| and [|T^. The data symbols 
are listed in the figure. 



1.47 GeV /c . The preliminary results are consistent with the Galster [|I^ parameterization 
from lower data 
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which, in parallel to the situation for Rp, implies much lower values of Geu in their Q range 
when coupled with Gmu values (either the precision data of [|1^] or the model fits). 

In this paper, in addition to the above comparison of GKex(Ol) and GKex(Ol-), we fit 
the model of type DR-GK'(l), with the added isoscalar vector meson a;'(1419) pole, to the 
following data sets, chosen to determine the effect of the old Geu and Gep data in direct 
conflict with the values of Rn and Rp from modern polarization measurements: 



The fit GKex(02L) from the full data set of [[l|] with the addition of [|T3| and [jT4|, the 
omission of [|| (as above for GKex(Ol-)) and the Geu values for > 0.779 GeV^/c^ 

of H, m, and m 



(b) The fit of GKex(02S) to the same data set as above except for the omission of the Gep 
values for Q^>1.75 GeVVc^ of §. 

It will be seen that the omission of the conflicting Geti data, GKex(02L), has a much bigger 
influence than the omission of GKex(Ol-), enabling a much better fit to Rp in addition 
to a very good fit to Rn, compared to GKex(Ol). With the removal of the conflicting Gep 
data, GKex(02S), the fit to all the remaining data, including Rp, is very satisfactory. 

In Section |I| we will specify the models and parameters used in this article, and the data 
sets used in Section |lll[ In Section |I^ we present the results of the four GKex fits in compar- 
ison with each other. We extrapolate beyond the present experimental range of momentum 
transfer where necessary for predicting available deuteron emff data. The model GKex(02S) 
fits the modern and consistent older data well and meets the requirements of dispersion 
relations and of QCD at low and high momentum transfer. Conclusions are presented in 
Section 0. 
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II. THE NUCLEON EMFF MODEL 

In fitting the nucleon emff data including the new Rn and Rp results we have chosen 
to use the extended GK model DR-GK'(l) of Ref. ^ with the addition of a pole term for 
the well established isoscalar vector meson Ci;'(1419), whose mass is lower than that of the 
already included isovector vector meson p'(1450). We denote this model as GKex. The choice 
of the particular parameterization DR-GK'(l) was made because of its low value and the 
fact that its predicted values of Rp were a little closer to the data than those of the other 
extended models. In addition DR-GK'(l) has the following good physical properties: 

(1) It uses the QCD cut-off A2 for the helicity flip meson- nucleon form factors, rather than 
the meson cut-off Ai used by DR-GK(3) and DR-GK'(3). 

(2) The evolution of the logarithmic dependence on is controlled by the quark-nucleon 
cut-off A^), along with Aqcd- DR-GK(l) and DR-GK(3) use A2 instead of A^- 

(3) Fitted to the data set of [||] it flnds Aqcd = -1163, close to the expected value. The 
form factors are not very sensitive to this parameter which is flxed at .15 for the flts 
to the new data sets. 

So that the reader need not make constant reference to we repeat the relevant formulas 
here together with the new ci;'(1419) terms. 

The emff of a nucleon are defined by the matrix elements of the electromagnetic current 

{N{p')\ J, \N{p)) = eu{p'){^,F,'^iQ') + ^a.^Q'^ F^^ {Q')]uip) (2) 

where N is the neutron, n, or proton, p, and —Q^ = {p' — pY is the square of the invariant 
momentum transfer. F^IQ"^) and F^IQ"^) are respectively the Dirac and Pauli form factors, 
normalized at = as 

Ff(0) = l, Fi"(0) = 0, F,'iO) = Kp, F^{0) = n^. (3) 

The Sachs form factors, most directly obtained from experiment, are then 

Gen(Q') = F.'^m-TFi^m 

GuM') = (Q^) + (Q^) . (4) 
Expressed in terms of the isoscalar and isovector electromagnetic currents 



2F[ = Fr + Fr, 2Fp = Fr-Fr, (^ = l,2) 



(5) 
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The GKex model has the following form for the four isotopic emff: 

_ , 1.0317 + 0.0875(1 + QVO.3176)-^ , 
(l + gVO.5496) ^ 

+ f^^2F'iQ") + (l - 1-1192 iV/2 - f)F^{Q^) 

_ ,, 5.7824 + 0.3907(1 + gVO.1422)-^ 

(l + gVO.5362) ) 

2 

+ Vt^^^T7W^2''(Q') + - 6.1731 iV/2 - Vt^)^2''(Q') 

Jp' '^p' + V ^ Jp' ^ 

+ (l-^-^)Ff(g^) (6) 

^ Juj Jul' ^ 

2 2 2 



. ( Qu Qijj' g(l,\ T?Dini'2 



Ju) Jl 

where the pole terms are those of the p, p', uo, uj' , and mesons, and the final term of 
each equation is determined by the asymptotic properties of PQCD. The F^, a = p, u, or (j) 
are the meson-nucleon form factors, while the Ff are effectively quark-nucleon form factors. 

For GKex the above hadronic form factors are parameterized in the following way: 



+g2VA2 + g2 



where a = p, and Ki^d is Ai for F", for F-^, 



withO^ - e,. to[(AUQ'^)/AgiCD] 
- <2 1„(AVA|„) ■ 

This parameterization, together with Eq. (^, guarantees that the normalization conditions 
of Eq. (Q) are met and that asymptotically 

i^i- [Q'ln(gVA|cD)]"' 

~ Fl/Q^ (8) 
i = is, iv 
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as required by PQCD. The form factor F{{Q'^) vanishes at = 0, and it and (Q^) 
decrease more rapidly at large than the other meson form factors. This conforms to the 
Zweig rule imposed by the ss structure of the meson [Q. 
This model has at most 14 free parameters: 

(i) 8 couplings to the pole terms, the 4 gm/ fm and the 4 for the p', u, u' and mesons. 

(ii) 4 cut-off masses in the hadronic form factors, Ai, A2, Ad, and /x^. 

iii) The mass determing the size of the logarithmic behavior, Aqcd- 

iv) The normalization factor N for the dispersion relation contribution of the p meson. 

However at most 12 of these parameters are freely varied in any of the fits made in the 
next section to the chosen data sets. 



III. DATA BASE AND FITTING PROCEDURE 

As previously stated, GKex(Ol) is the same as DR-GK'(l) of ||l[]. This model had the 
best fit to the full data set available at the publication of ||l|] with g'^/ = K''^ = ^ and with 
N=l. For GKex(Ol-) the 4 data points of were omitted from that data set. In this case 
dU fL and k'^ were still supressed but N was freely varied. 

In the fits GKex(02L) and GKex(02S) g^j/ and were freely varied, but these fits 
fixed N=l again (implying negligible error in the dispersion relation evaluation) and Aqcd 
was fixed at the physical value of 0.15 GeV/c. The important difference from the data 



set of GKex(Ol-) is the addition of the higher Q Rp data points of [0 and the i?„ data 
points of and the omission of the Geti values for > 0.779 GeV^/c^ of 0, |18[ and 
113. In the shorter data set of GKex(02S) the Gep values for > 1.75 GeVVc^ of § 
are also omitted. The free parameters were optimized using a Mathematica program that 
incorporates the Levenberg-Marquardt method. 



IV. RESULTS 

Table | presents the parameters which minimize the above 4 cases. For all 4 

parameter sets the hadronic form factor cut-off masses, Ai, A2, A^, and /x^. are reasonable. 
The relatively large value of A2, which controls the spin-flip suppression in QCD, is consistent 
with the slow approach to asymptopia observed in polarized hadron scattering. For the two 
cases in which Aqcd is a fitted parameter, as well as the two for which it is fixed, it is 
consistent with high energy experiment. The addition of the ci;'(1.419) meson in GKex(02L) 
and GKex(02S) has moved k^j closer to the expected small negative value than all earlier 
fits, but there is still the implication of some effect from a higher mass isoscalar meson. The 
adequacy of the fits is an indication that the form factors with more poles would be similar 
to those already obtained. 

In Table H the values of are listed for the 4 cases and the contribution from each of the 
six form factor classes of measurement are detailed. Also shown are the values of when 
any data points omitted from the fit are re-inserted. 

We note, as can also be seen in Figs. ^ and ^ that the quality of fit to the magnetic 
form factors, Gmp and Gmu changes negligibly as we refit to the datas sets that differ in the 
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TABLE I: Model parameters. Common to all models are = 3.706, Kg = —0.12, rup = 0.776 GeV, 
= 0.784 GeV, = 1.019 GeV, = 1.45 GeV and m^^ = 1.419 GeV. 



Parameters 


Models 




GKex(Ol) GKex(Ol-) 


GKex(02L) GKex(02S) 


a(o')/fto') 

il{h' J/ J (H J 


0.0636 


0.0598 


0.0608 


0.0401 




-0.4175 


-15.9227 


5.3038 


6.8190 


n If 

yojl Jul 


0.7918 


0.6981 


0.6896 


0.6739 




5.1109 


1.9333 


-2.8585 


0.8762 


n If 


— U.dUll 


— O.ozyU 


— U.i852 


— U.lbYb 




13.4385 


2.3241 


13.0037 


7.0172 




1.1915 


1.5113 


0.6848 


0.8544 








0.2346 


0.2552 








18.2284 


1.4916 


Ai 


0.9660 


1.1276 


0.9441 


0.9407 




1.3406 


1.8598 


1.2350 


1.2111 


As 


2.1382 


1.2255 


2.8268 


2.7891 


Aqcd 


0.1163 


0.1315 


0.150(^) 


0.150('') 


N 


i.o(^) 


0.8709 


l.O^*^) 


i.o(^) 



('^^not varied 



TABLE II: Contributions to the standard deviation, x^, from each data type for each of the 
models. The number of data points contributing is in parentheses. For each data type the first row 
corresponds to the data set for which the model parameters were optimized, the second row to the 
full data set. 



Data 


Data 




Models 






type 


set 


GKex(Ol) 


GKex(Ol-) GKex(02L) GKex(02S) 


Gmp 


opt 


43.3(68) 


43.6(68) 


48.1(68) 


47.9(68) 




full 




same as above 






Gep 


opt 


67.2(48) 


48.2(44) 


75.3(44) 


30.5(36) 




full 


67.2(48) 


74.8(48) 


112.2(48) 


136.8(48) 




opt 


122.4(35) 


120.2(35) 


121.0(35) 


122.7(35) 




full 




same as above 






Geu 


opt 


64.8(23) 


64.2(23) 


24.1(15) 


24.2(15) 




full 


65.3(24) 


65.0(24) 


68.2(24) 


68.3(24) 


Rp 


opt 


29.0(17) 


22.6(17) 


23.1(21) 


11.8(21)1 




full 


114.0(21) 


106.5(21) 


23.1(21) 


11.8(21) 


Rn 


opt 


0.0(0) 


0.0(0) 


0.6(3) 


0.6(3) 




full 


9.6(3) 


17.7(3) 


0.6(3) 


0.6(3) 


Total 


opt 


326.7(191) 


298.9(187) 


336.3(195) 


237.7(178) 




full 


421.8(199) 


427.8(199) 


369.2(199) 


388.1 
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FIG. 3: Gmp normalized to UpGd- Comparison of the models GKex(Ol) [solid], GKex(02L) [dotted] 
and GKex(02S ) [dash-dotted] with the data of [|, [§, [|], [|], [0, |l|], [|l| and [|^. (a) The 
full data range, (b) Expansion of the range < 3.0 GeV'^/c'^. The data symbols are listed in the 
figure. 



electric form factors and the electric to magnetic form factor ratios. As discussed in [|T|, the 
large excess of over the number of data points for Gmu is due to obvious inconsistencies 
in the data set for Gmu at < 0.8 GeV/c^. The displacement of adjacent data points 
well beyond their error bars in this range is evident in the figures and contributes about 90 
to the of Gmti- 

The interesting changes are, of course, in the fits to Gep, Geu, Rp and As noted in 
the introduction, removing the 4 very high values of Gep data [^ does surprisingly little to 
allow a better fit to the Rp data already in GKex(Ol). Several of the parameters, all three 
Krai A2 and Kdi have large changes (see Table |), but this results in a small shift between 
the predictions of GKex(01 and GKex(Ol-) as is evident in Table |T| and Figs. ^ and 0. The 
figures show a slightly better fit to Rp correlated with a very slightly worse fit to Gep- The 
former is reflected in Table |T| by the decrease in the contribution of the 17 Rp points 
to which those cases were optimized from 29.0 tO 22.6. When the 4 higher of [jTSf are 
added the contribution is much larger than the number of points (21). The drop in the 
X^ contribution to Gep from 67.2 to 48.2 is entirely due to the omission of the 4 data points 
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FIG. 4: Gmu normalized to fJ-nGd- Comparison of the models GKex(Ol) [solid], GKex(02L) [dotted] 
and GKex(02S ) [dash-dotted] with the data of [|, |l|], Q, [|l8|, ||, |2|, [||, [|4|, |||, |2§, 
and |2^. The data symbols are listed in the figure, (a) The full data range, (b) Expansion of the 
range Q'^ <2.0 GeV^/c^. 



of [^, but the for the full set of 48 points is a little larger because of the competition with 
Rp. The implication is that there is a constraint on the fit to Rp from data independent of 
Gep, arising from the model correlations between all the nucleon emff. This is shown to be 
the case below. 

Substituting the new Rn values for the conflicting Geu data of 
a large difference between the GKex(02L) and GKex(Ol-) fits to Gep, Ge 
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.p and Rn 

seen in Table || and Fig. § - Fig. p. In particular for GKex(02L) the contribution for all 
21 Rp data points is 23.1. Fig. ^ shows the strong improvement in the fit to Rp. The figure 
also shows that the goodness of the value is somewhat misleading because that fit is 
systematically high for the 3 highest data points. On the other hand the contribution 
for all 44 Gep data points increases from 48.2 in GKex(Ol-) to 75.3 in GKex(02L) because of 
the compromise of better fitting Rp. The contribution for the 3 Rn points now included 
is only 0.6. Geu now contributes 24.1 for the remaining 15 data points (which still include 
highly scattered low data as discussed in []I|) instead of 64.8 for the 23 data points in 
GKex(Ol-). 

For the GKex(02S) case the Gep data of which is clearly inconsistent with the new 
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is also omitted. The results are very good if the modern data is chosen 
when in conflict with the older Rosenbluth separation results. The contribution from the 
remaining 36 Gep points is only 30.6 and for all 21 Rp points onlyll.8. For the remaining 
types of form factor measurements there is a negligible change of between the GKex(02L) 
and GKex(02S) cases. 
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FIG. 5: Gep normalized to Gd- Comparison of the models GKex(Ol) [solid], GKex(02L) [dotted] 
and GKex(02S ) [dash-dotted] with the data. The data references are the same as in Fig. 1 and the 
data symbols are listed in the figure. The points labelled by open circles are obtained by multiplying 
Rp data Q and [|l^ by the Gmp of GKex(02S) normalized by fipG^- (a) The full data range, (b) 
Expansion of the range < 2.0 GeV^/c^. 



Figs. 1^, 0, and || show the successive improvements in Gep, Rp, Geu, and i?„ as 
the optimization data sets are varied from GKex(Ol) to GKex(02L) and to GKex(02S). To 
demonsrate the correlation between the electric form factors and the ratio of electric to mag- 
netic form factors we have, in Figs. |^ and |^, entered (as circles) the electric form factor values 
obtained by multiplying the experimental Rp and Rn values by the case GKex(02S) model 
values of the magnetic form factors normalized by the magnetic moments. The correlation 
with the model prediction for the electric form factors is excellent. 

The figures show the model extrapolations of Rp, Gmti, Geu and Rn up to of 8 GeV^/ 
for the guidance of future experiments and because of their relevance to deuteron and other 
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FIG. 6: Rp, the ratio fipGEp/GMp- Comparison of the models GKex(Ol) [sohd], GKex(02L) [dotted] 
and GKex(02S ) [dash-dotted] with the data. The data is from |5|, and p^ . 

nuclear electromagnetic scattering predictions. The extrapolation is sensitive to the weight 
given to the polarized vs. the Rosenbluth separation data in the fits. The resolution of this 
dichotomy will, in the context of the physical model employed here, greatly restrict the 
nucleon emff over a large range of momentum transfer. 

The polarization measurements of Rp and Rn may soon be extended to larger Q^, so 
it is of interest to examine the predictions of the good fit GKex(02S) as increases. 
As seen in Fig. ^ the model curve is, as is the data, approximately linear in the range 
0.5 GeV^/c^ < < 5.6 GeV^/c^, but the model curve's slope is gradually decreasing in 
manitude. A linear fit to the data would change sign near = 8 GeV^/c^ where the model 
predicts 0.19. The model crosses zero near = 14 GeV^/c^ with a very small slope. 

Fig. ^ shows the Galster curve, -R^"'**'^'"((5^) of Eq. (|l|), to compare with the model and 
the data. The model fits the data but deviates from the Galster curve after that. The model 
increases faster, reaching a maximum of 0.426 at Q2 = 4 QeV^/c^ where the Galster value 
is only 0.309, while increases monotonically to an asymptotic value of 0.342. A 

measurement of the present quality at = 4 GeV^/c^ could distinguish between the two. 

V. CONCLUSIONS 

The GKex model, consistent with vector meson dominance and perturbative QCD in 
the appropriate momentum transfer regions, represents well a consistent set of neutron and 
proton emff. This set includes polarization measurements, which are directly related to the 
ratios of electric to magnetic form factors, and differential cross section measurements of the 
magnetic form factors. The values of the electric form factors from the Rosenbluth separation 
of the differential cross section is, in our final selection GKex(02S), only used for the lower 
range of where the magnetic contributions are less dominant. Because of the physical 
properties of the model and the good quality of the fit we expect that the model predictions 
are sufficiently accurate to be used for predictions of the electromagnetic properties of nuclei. 
The model values may also be useful in planning future experiments. 

The above conclusions are only valid to the extent that adequate physics is included 
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FIG. 7: normalized to Gd- Comparison of the models GKex(Ol) [solid], GKex(02L) [dotted] 
and GKex(02S ) [dash-dotted] with the data of [|, 0, (l|], [|§, |§, |2|, [|g, [0, ||, g, 
||, H and [||. The data of [|o|, [||] and [||] are the reevaluated values of The slope at 
= is from [^. The points labelled by open circles are obtained by multiplying i?„ data [14] 



by the Gmu of GKex(02S) normalized by fJ-nGd- (a) < 8.0 GeV^/c^. (b) Expansion of the range 
Q2 < 2.0 GeVVc^. 



in the GKex models. Only the p meson exchange includes the width of the vector mesons 
(from dispersion relation results). There will be corrections from the widths of the other 
exchanged vector mesons. However the next most important, the u and </>, have very narrow 
widths. The higher masses of the p' (1450) and the u' (1420) reduces the importance of 
their substantial widths because of their distance from the physical region and their partial 
replacement by the pQCD term. 

In assuming vector dominance we have neglected the multi-meson exchange continuum 
contributions. In particular the two-pion continuum may have an influence at very low 

< 0.4 GeV^/c^. Indeed, as remarked in[|T| and can be seen in Fig. 5b, the Gep data of [^ 
has a more negative slope for < 0.3 GeV^/c^ than the higher data and the model fit. 
The addition of a 2-pion exhange term to the model may enable a change of slope between 
the two regions, but would have little effect on the model fit for > 0.5 GeV^/c^. 

One may also want to consider some higher mass vector mesons. This would have some 
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* [14] 



FIG. 8: Rn, the ratio UnGEn/GMn- Comparison of the models GKex(Ol) [solid], GKex(02L) [dotted] 
and GKex(02S ) [dash-dotted] with the data. The dashed curve is R^''^'^'"' (Q^) of Eq. @. The 
"experimental" points are described in the text [p^] . 

importance in the fits of [Q] and|^, but are much less important in these GK type models 
because of the transition to pQCD behavior. 
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